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(54) Suppression of undesired components in measured spectra 



(57) The effects of unwanted components such as 
H 2 0 and C0 2 in spectral data measured by a spectrom- 
eter such as an FT-IR spectrometer are suppressed by 
a technique in which data representing high resolution 
spectra of the unwanted component is acquired, this 
data is modified so that its resolution matches that of the 
instrument, the results data is filtered to allow for per- 
turbing effects of the sample and the resulting data is 
subtracted from the measured sample spectrum to pro- 
vide corrected output data. 
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Description 

[0001] This invention relates to spectrometers and in 
particular relates to suppressing undesired spectral 
components in the spectra obtained from the spectrom- 5 
eters. The invention has particularly but not exclusive 
application to FT-IR spectroscopy. 
[0002] In for example an FT-IR spectrometer, infra-red 
or near infra-red radiation is directed from a source of 
such radiation towards a sample under investigation. 10 
Radiation transmitted by or reflected from the sample is 
received at a detector or receiver and the output of the 
detector is processed by a signal processor in order to 
obtain the spectral characteristics of the sample. In car- 
rying out measurements it is first necessary to obtain is 
what is known as a background measurement that is to 
say to the measure the background spectrum without a 
sample in place at the sample station. Subsequently 
measurements are made with the sample in place and 
the desired sample spectrum is obtained from the ratio 20 
of the measurement obtained with the sample in place 
to the background measurement. 
[0003] A significant proportion of spectroscopy carried 
out in the mid infra-red range is concerned with spectra 
of solids and liquids which have absorbtion bandwidths 25 
in the tens of wave numbers. As a consequence, many 
of the measurements are made at, for example, 4cm-1 
or 2cm-1 resolution. At such moderate resolutions water 
vapour having line widths closer to 0.1cm-1 is strongly 
under resolved and bands that might otherwise be satu- 30 
rated with absorbance well in excess of 1.0 are broad- 
ened to the point where they show relatively low peak 
absorbance. Two major consequences result from this 
under resolution: The absorbance becomes substan- 
tially non-linear with concentration causing the absorb- 35 
ance spectrum shape to become a strong function of 
concentration and in addition the lineshapes become 
entirely dominated by the instrument lineshape function 
making the spectra have a different character according 
to the instrument type and set up. More significantly 40 
since the instrument lineshape function may well be 
influenced by the sample or sampling accessory 
through vignetting another beamed geometry distur- 
bances the water vapour spectrum in the sample spec- 
trum may not entirely resemble the nomunally similar 45 
spectrum in the background. That is to say the effects 
do not cancel out when the sample spectra and the 
background spectra are ratioed. The overall result is 
that it is extremely difficult to subtract out the effects of 
water vapour consistently by any means of linear spec- so 
trum differencing typically employed in this type of spec- 
troscopy even if advance statistical methods such as 
Principal Components Analysis are employed. A similar 
problem exists with other unwanted components such 
as those of carbon dioxide. 55 
[0004] A further aspect of the problem is related to the 
overlap of features in the sample spectrum with features 
in the water vapour spectrum. Even when an exact 



spectrum of water vapour measured under the current 
sampling conditions can be generated independently by 
some means, it is very difficult to estimate precisely the 
proportion of water vapour spectrum that must be sub- 
tracted out from the detected spectra especially by 
automatic algorithmic approaches. 
[0005] The present invention is concerned with a 
novel technique for subtracting out from the measured 
sample spectrum undesired components such as those 
arising from water vapour and carbon dioxide. 
[0006] According to the present invention there is pro- 
vided a spectrometer which comprises a source of ana- 
lysing radiation, a detector for detecting radiation 
transmitted through or reflected from a sample under 
investigation, and a processing means for processing 
the output of the detector to produce spectral data relat- 
ing to a sample under investigation, wherein the proces- 
sor is arranged, in order to suppress the effects of 
unrequired components such as water vapour or carbon 
dioxide in the spectral data, to acquire reference data 
representing high resolution spectra of the unwanted 
component or components, to modify said data so that 
its resolution simulates that of the spectrometer, to filter 
said data so as to allow for pertutbing effects in the sam- 
ple spectrum and to subtract the resulting data from the 
measured sample spectrum in order to provide a cor- 
rected output data. The acquired data may represent 
said high resolution spectra at a plurality of tempera- 
tures. The modification may comprise processing said 
reference data so as to broaden the resolution to that 
which matches the sample spectral data being meas- 
ured. The broadening may be carried out by a convolu- 
tion technique in which the reference data are 
convolved with a computed line shape function appro- 
priate to the spectrometer. These steps may be 
repeated for different reference data representing differ- 
ent concentrations of unwanted component or compo- 
nents. 

[0007] The processor may also be arranged to gener- 
ate perturbed versions of the broadened reference data 
to take into account variations in at least some operat- 
ing parameters such as temperature and optical line 
width. 

[0008] The filtering may include creating a filter for fil- 
tering the spectra to emphasise higher resolution parts 
of the unrequired components. The filter may comprise 
a band pass filter. Filtering may also include carrying 
out a least square fit of the filtered sample spectral data 
to the filtered unwanted component spectra. Addition- 
ally the processor may iterate the least square fit to 
remove those parts of the spectrum which have a poor 
fit. 

[0009] The processor may be arranged to use the 
resulting co-efficients to compute an unfiltered 
unwanted component spectrum. This is then subtracted 
from the measured sample spectrum in order to obtain 
a corrected sample spectrum. 
[0010] The invention will now be described now by 



3 



3 



EP 0 982 582 A1 



4 



way of example and with particular reference to the 
accompanying drawings. In the drawings: 

Figure 1 is a schematic illustration of an infra-red 
spectrometer in which the present technique can be 5 
implemented and Figure 2 is a flow chart illustrating 
the sequence of steps carried out by a processor 
arranged to operate in accordance with the present 
invention. 

10 

[0011] Referring to Figure 1 of the drawings an infra- 
red spectrometer of the single beam type such as an 
FT-IR spectrometer in very general terms, comprises a 
source of analysing radiation 10 which is arranged to 
direct analysing radiation to a sample station 12. Typi- is 
cally the source of analysing radiation will be associated 
with an interferometer of the Michelson type which pro- 
duces a scanning beam by means of which the sample 
can be irradiated. Radiation from the sample station 12 
is received at a detector or receiver 14 and the output of 20 
the receiver is processed by a processor 16 to provide 
data representative or the spectrum of a sample under 
investigation. Generally speaking a spectrometer of this 
type has an associated PC for carrying out at least 
some of the processing and enabling an operator to 25 
issue appropriate instructions to operate the spectrom- 
eter. The way in which a spectrometer of the type shown 
in Figure 1 operates in order to obtain a spectrum will be . 
known to those skilled in the art and will therefore not be 
described in any more detail since it is unnecessary for 30 
understanding the present invention. 
[0012] Embodiments in accordance with the present 
invention are concerned with routines which are oper- 
ated by the processor in order to suppress or cancel out 
the effects of undesired components such as water 35 
vapour or carbon dioxide in the measured sample spec- 
tra. It is envisaged that such routines can be operated 
by a wide range of spectrometers provided that the 
processing means of the spectrometer has sufficient 
capability to run the software involved in carrying out the 40 
routines. The description given here is directed in the 
main to the way in which the invention can be imple- 
mented in FT-IR spectroscopy. 
[001 3] Initially a general description of the method will 
be given followed by an implementation in more specific 45 
terms. The description will refer to removing the effects 
of water vapour but it should be understood that the 
technique applies equally well to supparesion of other 
undesired components such as carbon dioxide. 
[0014] In very general terms the processor is arranged so 
to remove or suppress the effects of any water vapour 
components in the spectrum by the following steps. 

1. Generating a properly broadened approximately 
representative set of spectra for the unwanted com- 55 
ponent which span the range of non-linear behav- 
iour expected in the instrument. 

2. Generating perturbed spectra to take into 



account a degree of variability in lineshape param- 
eters; and 

3. Iterative fitting using a combination of filtering 
techniques. 

The initial step in the method is to acquire stored high 
resolution spectral data relating to water vapour. This 
should have sufficient resolution to accurately charac- 
terize the natural lineshape. It is only when the line 
shape is properly resolved that the behaviour of the 
bands can be properly represented as peak absorb- 
ances increase and the non-linearity of the transmission 
values becomes significant. Such data is available com- 
mercially and one example which can be used is a set 
of peak tables known as (HITRAN). Significant parame- 
ters of such spectral data are atmospheric pressure 
which under normal circumstances is responsible for 
line broadening and in the present case will be under- 
stood to be one atmosphere, and a representative gas 
temperature which in the present case is taken to be 
23°C, a typical laboratory environment. 
[0015] Starting from this single spectrum the method 
generates spectra for three different concentrations of 
water at a pathlength typical of the actual pathlength 
encountered in practice eg 25cm. The concentrations 
selected are 25% relative humidity, (RH) 50%RH and 
100%RH thereby spanning the range of concentrations 
typically encountered in practice. The simple way of 
obtaining these spectra is to square the transmittance 
spectrum of a base 25%RH spectrum to obtain the 
50%RH spectrum and then repeat the process in order 
to obtain the 100%RH spectrum. An alternative way of 
doing this would be by a conventional logarithmic scal- 
ing technique. 

[0016] These three principal spectra are then broad- 
ened in transmittances in order to simulate approxi- 
mately the broadening that occurs in the spectrometer 
itself. In the present case, ie a Fourier Transform spec- 
trometer, the spectrum is broadened on a linear wave 
number abscissa using an FT lineshape function arising 
from a combination of the interferogram scan length 
used in the instrument and an apodisation function 
employed. The broadening can be performed either in 
spectrum space with a convolution or interferogram 
space by multiplying an envelope function. In the 
present example the operation is carried out in the 
spectrum space. In the next step the spectrum is broad- 
ened on a logarithmic wave number abscissa using a 
rectangular line shape function associated with a per- 
fectly behaved circular Jacquinot stop. The Jacquinot 
stop controls the divergence of rays passing through the 
interferometer of the spectrometer and effectively con- 
trols its resolving power. The associated line width is 
proportional to wave number and hence on a logarith- 
mic scale appears to have constant width. There are a 
number of ways in which this calculation can be per- 
formed. In the present case it is achieved by computing 
the mean value over the line width from the area under 
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spectrum, using a variable line width. It should be 
understood that the Jacquinot stop can also introduce a 
spectrum shift. In the present case this is already sub- 
stantially calibrated out by use of a centre broadening 
function. 

[0017] At this stage the process has produced three 
principal broadened spectra at the various concentra- 
tions. The next step is to generate three additional minor 
components which allow for variation from the com- 
puted spectra. The first is a similarly broadened spec- 
trum at 50%RH but at 10°C higher temperature to allow 
for temperature variation. The second is at 50%RH at 
the standard temperature broadened with a slightly 
smaller Jacquinot stop function to allow for imperfect 
optics. The third is a spectrum shifted in wave number 
by a small constant to allow for calibration errors. In 
each case, the minor component actually employed in 
the subsequent processing is the difference between 
the perturbed spectrum and the principal 50%RH com- 
ponent. 

[0018] In the majority of cases the six components 
that it is to say the three major components and the 
three minor components should be sufficient to imple- 
ment a technique in accordance with the present inven- 
tion. However, if necessary it would be possible to 
iterate the fit using more appropriate starting parame- 
ters once estimates for the parameters have been made 
from the fit components, for example temperature and 
path length. In some cases it may be possible to operate 
with less than six components especially if an iterative 
approach of refining the initial parameter estimates is 
adopted. An iterative approach works best for single 
beam spectra where the water vapour spectrum is 
directly represented as opposed to a ratio spectrum 
where only small differences in water vapour may be 
present. Once having refined the parameters from the 
single beam subsequent elimination of the effects of 
water vapour is best carried out on the ratio spectrum 
since the ration itself eliminates a significant part of the 
water vapour effect. 

[001 9] Having established the components which can 
be used to accurately reconstruct the expected water 
vapour trans mitta nee spectrum, it is necessary to deter- 
mine the proper coefficients for the fit in the presence of 
the unknown and potentially significant determining fac- 
tors in the data arising from, for example, the spectrum 
of the sample. These may be many times stronger than 
the spectrum of the water vapour. It is found in practice 
that the problem of subtracting out the water vapour 
remains linear in absorbence even though the construc- 
tion of the water vapour spectrum itself is a non-linear 
problem. Thus, in the present example the spectra are 
converted to absorbence before the following opera- 
tions. 

[0020] The first step in reducing the interference from 
the sample spectra is to filter all the data, that is the cal- 
culated component spectra and the sample spectrum, 
with a band pass filter that accepts data centred around 



10cm-1 resolution while rejecting data at higher and 
particularly lower resolution. The parameters of this fil- 
ter are based on the general properties of infra-red 
spectra of solids and liquids typically measured in the 
5 applications of concern. The underlying principle is that 
of a matched filter. 

[0021] In order to create such a filter typical sample 
spectra obtained from a commercial library are Fourier 
transformed and the square root sum of squares of the 

10 real and imaginary components calculated in order to 
give an idea of the distribution of the information present 
across the range of resolution. In the present case this 
indicates the amount of data at resolutions better than 
10cm-1 is relatively small. A similar calculation for water 

15 vapour shows that being a gas there is considerable 
information up to the resolution limit of the spectrum. 
The ratio of the data for the water vapour versus the 
library data shows a broad peak centred around 10cm- 
1. The shape of this peak approximates to a gaussian 

20 and this function forms the basis of the matched filter, 
and is analytic in both transform spaces. 
[0022] It should be noted that the matching of the filter 
is not especially critical and good results can be 
obtained using a simple first derivative filter. However, a 

25 gaussian function approximates closely to the ideal for a 
general case. 

[0023] The above described resolution filtering 
reduces the effects of interference to water vapour least 
square fits significantly. However, it is still not sufficient 
30 to produce a satisfactory result in many situations. In 
order to achieve this an additional step of iterative 
weighting of the fit in the spectrum domain is carried 
out. 

[0024] In this the initial weight function employed is the 

35 measured transmittance spectrum itself. The reason for 
this is that regions of low transmittance have poor signal 
to noise ratio and a good fit cannot be expected in these 
regions. The weighted least squares fit of the filtered 
data is computed and the residual spectrum calculated. 

40 The weighted root mean square residual value is calcu- 
lated and then the weight function is modified to set zero 
weight for those portions of the spectrum where the res- 
idue exceeds four times the rms value. This process 
effectively blanks out regions of the spectrum where the 

45 fit is particularly bad. The fit is iterated and the process 
repeated until the fit is sufficiently improved. There may 
be a compromise between the number of iterations and 
the severity of the blanking and in practice it is found 
that the final result is better with less severe blanking 

so and more iterations. 

[0025] Having calculated the best parameters of the fit 
using filtered data the process reconstructs the full 
water vapour spectrum using the same fit parameters 
with the unfiltered data. The measured sample spec- 

55 trum is then divided in transmittance by the water 
vapour spectra in order to elimininate the water vapour 
component. 

[0026] It is believed that the above description will 



5 



7 



EP 0 982 582 A1 



8 



enable a reasonably competent computer programmer 
to implement a program for carrying out the process 
described. 

[0027] In order to assist in this respect the following 
description in conjunction with the flow chart of Figure 2 5 
is included as a specific example for eliminating an 
unwanted component in accordance with the present 
invention. 

[0028] This program effectively has two inputs these 
being as follows: 10 

1 The first is the target spectrum or the measured 
spectrum of the sample from which the water 
vapour component is to be subtracted. The spec- 
trum is typically a low resolution spectrum meas- 15 
ured at around 4cm-1 resolution in transmittance, 
digitised at a constant wave number interval at 
around 2x over-sampling, that is at 1cm-1 interval. 

2 This is the high resolution reference spectra of 
water vapour at two difference temperatures. These 20 
spectra typically resolve the natural line width of the 
water vapour under measurement conditions of the 
target spectrum, ie around 0.1cm-1 for normal pres- 
sure, and are stored in transmittance at intervals of 
1/64cm-1. The two temperatures should lie on 25 
either side of the temperature of the target spec- 
trum and the water vapour concentration should be 

set around 25%RH at the path length of the target 
spectrum. 

30 

[0029] Referring now to Figure 2 the first step illus- 
trated at 20 is the reading of the target spectrum that is 
to say the spectrum measured by the spectrometer. The 
next step illustrated at 21 is to read in the lower temper- 
ature reference spectrum, ie the second input referred 35 
to above. 

[0030] The next step shown at block 22 is to compute 
the lineshape corresponding to the Fourier transform by 
which the target spectrum was generated. This line 
shape is calculated by computing either analytically or 40 
numerically the Fourier transform of the apodisation 
function applied to the interferogram from which the tar- 
get spectrum was generated. In the case of no apodisa- 
tion the lineshape will be a sine function of an 
appropriate width. 45 
[0031] The next step shown at block 23 is to broaden 
the reference spectrum to the same nominal resolution 
as the target spectrum by convolving with the computed 
lineshape. Convolving involves multiplying the first part 
of the reference spectrum by a short line function and so 
summing the products. Then the alignment of the two 
spectra is shifted along by one point and the process 
repeated. The sums form a new spectrum which is the 
convolution of the original spectrum with the lineshape 
function. 55 
[0032] The next step is to reduce the sampling of the 
broaden reference spectrum to the same interval as the 
target spectrum by dropping data points. It should be 



noted that once the lineshape has been broadened data 
points can be discarded without loss of information. 
[0033] The next step, as illustrated at 25, is to square 
the original low temperature reference spectrum to pro- 
duce a 50%RH spectrum and this is broadened in a 
manner similar to that described for the initial reference 
data. 

[0034] Referring to block 26 this is then repeated to 
give the 100%RH broadened reference spectrum. At 
this stage there exist three spectra for that cover the 
concentration range of the target spectrum. The 
50%RH spectrum will be referred to as the principal 
component. 

[0035] The next step as illustrated at 27 is to read in 
the higher temperature reference spectrum. This is 
squared and broadened in a manner similar to that 
already described to a give a broadened 50%RH higher 
temperature fourth component. This will be used later to 
allow for any temperature variation in the target spec- 
trum. Each of the four components are then further 
broadened by a rectangular shape function which repre- 
sents the ideal effect of the spectrometers optical aper- 
ture. This rectangular function has a theoretical 
fractional optical frequency width 

dnu/nu = l-cos(theta) where that is the semi-angle 
subtended by the Jacquinot stop aperture at the colli- 
mating mirror. 

Theta=arctan (j-stop diameter/2* focal length) 

[0036] The convolution is computed directly by numer- 
ical integration of the area under this rectangle when 
multiplied into the operand spectrum. 
[0037] At this stage the process has computed four of 
the six required component spectra. The next step is to 
generate in a similar manner a j-stop broadened version 
of the principal component but with 1 0% narrower rec- 
tangular functions. This allows for some experimental 
variability of optical line width. 
[0038] The principal component is then interpolated to 
generate a sixth component with a small wave number 
shift proportional to optical frequency. This allows for 
some tolerance in the frequency calibration of the tar- 
get. 

[0039] Then all six components are converted to 
absorbance as illustrated at block 30 and the principal 
component subtracted from the last three so that the 
process handles only minor perturbations produced by 
temperature, j-stop and shift. 

[0040] The next step, shown in 31 is to produce a 
guassian-shaped band pass digital filter function of 
width (1/e) of 0.25 the Nyquist frequency and centred at 
0.25 the Nyquist frequency. This filter can be generated 
analytically and serves to emphasize high resolution 
information unique to the water vapour spectrum. In the 
next step 32 this filter is convolved with all components 
and also with the target spectrum. Then a gaussian filter 
is computed which has a similar width which is centred 
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at zero frequency (33). This filter is convolved with the 
unfiltered principal component to generate a weight 
spectrum (34) The filter so generated effectively 
smooths the absorbance peaks of the water vapour 
spectrum to generate a function that resembles the 
envelope of the peaks. The weight function is further 
multiplied by the unfiltered transmittance target spec- 
trum. 

[0041] The next step illustrated at 35 is to carry out a 
least squares fit of the six filtered components to the fil- 
tered target data using the weights. This composite 
weight function ensures closeness of fit wherever the 
water vapour is likely to be most intrusive. 
[0042] The next step shown by block 36 is to compute 
the square of the residual spectrum multiplied by the 
weight spectrum. The residual spectrum is the differ- 
ence between the target and the fitted spectrum. The fit- 
ted spectrum is computed by summing the component 
spectra scaled by the coefficients derived in the least 
squares fitting process. 

[0043] Referring to block 37 the next step is to com- 
pute the root weighted mean square residue. This 
number is the analog of the rms residue in an 
unweighted fit. 

[0044] Then the process computes a weighted resid- 
ual spectrum as illustrated at block 38 and then identi- 
fies data points in the weighted residual spectrum which 
are greater than 4x the root weighted mean square res- 
idue and zeros the corresponding weights in the weight- 
ing spectrum. This process identifies data points in the 
target that are beyond the acceptable range and dis- 
cards them as illustrated at 39. 
[0045] The process is then iterated (40) from the least 
square fit step typically six times in order to refine the fit 
and to make sure that the majority of the unacceptable 
data has been removed. 

[0046] The process then uses the fit coefficients gen- 
erated in the final iteration to compute an equivalent of 
the fitted spectrum from the unweighted unfiltered com- 
ponents. (See block 41.) 

[0047] This reconstructed water vapour spectrum is 
then subtracted from the unfiltered target spectrum (42). 
The reconstructed spectrum should accurately match 
the water vapour component in the target spectrum 
which should then subtract out cleanly. The result is the 
sample spectrum without the water vapour component. 
[0048] By way of general comments the effects of 
water vapour are most prominent in a well separated 
spectrum regions for example the regions centred 
around 1600cm-1 and the region around 3700cm-1. In 
theory since these are ail part of the same water vapour 
spectrum, one set of fit parameters should suffice for all 
regions. However modelling of these spectrum may not 
be exact and it may be necessary to treat the regions 
independently. 

[0049] At higher optical frequencies the instrument 
line shape function can become increasingly dominated 
by geometric and optical factors in the beam path. The 



method described above relies for its basis spectra on a 
reasonably close approximation to ideal behaviour 
since the Jacquinot stop is assumed to give rise to a 
rectangular line shape. Although the fit is not particu- 

5 larly sensitive to exact details of this line shape, being 
dominated for the most part by FT line shape, it may be 
preferable in some situations to use an alternative 
broadening function for the Jacquinot stop effect. 
[0050] The parameters of the least squares fit can 

10 yield diagnostic information on instrument performance. 
For example the overall strength of the water bands in 
the background spectrum can be used to indicate the 
state of purge of the instrument, the amount of shifted 
spectrum employed in the fit can be used as a check on 

15 the abscissa calibration. It could in principle be used as 
the primary abscissa calibration, although this may not 
be wholly desirable. 

[0051] The method described above relates to sup- 
pression of water vapour components. The principles 

20 described can be applied equally well to removal of car- 
bon dioxide components. However, account needs to be 
taken of differences arising from the difference in spec- 
trum and resolution distribution. Carbon dioxide has a 
more regular band strength structure that is easily 

25 obscured by underresolving but on the other hand it has 
less overlap with typical sample spectra, particularly in 
the strong 2350cm-1 band. This tends to place more 
emphasis on the spectrum weighting and less on reso- 
lution filtering. 

30 

Claims 

1. A spectrometer which comprises a source of ana- 
lysing radiation, a detector for detecting radiation 

35 transmitted through or reflected from a sample 
under investigation, and a processing means for 
processing the output of the detector to produce 
spectral data relating to a sample under investiga- 
tion, wherein the processor is arranged, in order to 

40 suppress the effects of unrequired components 
such as water vapour or carbon dioxide in the spec- 
tral data, to acquire reference data representing 
high resolution spectra of the unwanted component 
or components, to modify said data so that its reso- 

45 lution simulates that of the spectrometer, to filter 
said data so as to allow for perturbing effects in the 
sample spectrum and to subtract the resulting data 
from the measured sample spectrum in order to 
provide a corrected output data. 

50 

2. A spectrometer according to claim 1 wherein the 
acquired data represent said high resolution spec- 
tra at a plurality of temperatures. 

55 3. A spectrometer according to claim 1 or claim 2 
wherein the modification comprises processing 
said reference data so as to broaden the resolution 
to that which matches the sample spectral data 
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being measured. 

4. A spectrometer according to claim 3 wherein the 
broadening is carried out by a convolution tech- 
nique in which the reference data are convolved 5 
with a computed line shape function appropriate to 
the spectrometer. 

5. A spectrometer according to claim 4 wherein these 
said steps are repeated for different reference data 10 
representing different concentrations of unwanted 
component or components. 

6. A spectrometer according to any one of claims 3 to 

5 wherein the processor is arranged to generate 15 
perturbed versions of the broadened reference data 
to take into account variations in at least some 
operating parameters such as temperature and 
optical line width. 

20 

7. A spectrometer according to any preceding claim 
wherein the filtering includes creating a filter for fil- 
tering the spectra to emphasise higher resolution 
parts of the unrequired components. 

25 

8. A spectrometer according to claim 7 wherein the fil- 
ter comprises a band pass filter. 

9. A spectrometer according to claim 8 wherein the fil- 
tering also includes carrying out a least square fit of 30 
the filtered sample spectral data to the filtered 
unwanted component spectra. 

10. A spectrometer according to claim 9 wherein the 
processor is also arranged to iterate the least 35 
square fit to remove those parts of the spectrum 
which have a poor fit. 

11. A spectrometer according to claim 9 or claim 10 
wherein the processor is arranged to use the result- 40 
ing co-efficients to compute an unfiltered unwanted 
component spectrum. 

12. A method of operating a spectrometer which com- 
prises a source of analysing radiation, a detector for 45 
detecting radiation transmitted through or reflected 
from a sample under investigation, and a process- 
ing means for processing the output of the detector 

to produce spectral data relating to a sample under 
investigation, wherein method being carried out in so 
order to suppress the effects of unrequired compo- 
nents such as water vapour or carbon dioxide in the 
spectral data, and comprising the steps of acquiring 
reference data representing high resolution spectra 
of the unwanted component or components, to 55 
modifying said data so that its resolution simulates 
that of the spectrometer, to filtering said data so as 
to allow for perturbing effects in the sample spec- 



trum and to subtracting the resulting data from the 
measured sample spectrum in order to provide a 
corrected output data. 

13. A method according to claim 12 wherein the 
acquired data represent said high resolution spec- 
tra at a plurality of temperatures. 

14. A method according to claim 12 or claim 1 3 wherein 
the modification step comprising processing said 
reference data so as to broaden the resolution to 
that which matches the sample spectral data being 
measured. 

15. A method according to claim 14 wherein the broad- 
ening is carried out by a convolution technique in 
which the reference data are convolved with a com- 
puted line shape function appropriate to the spec- 
trometer. 

16. A method according to claim 15 including repeating 
said steps for different reference data representing 
different concentrations of unwanted component or 
components. 

17. A method according to any one of claims 12 to 16 
including generating perturbed versions of the 
broadened reference data to take into account vari- 
ations in at least some operating parameters such 
as temperature and optical line width. 

18. A method according to any one of claims 12 to 17 
wherein the filtering step includes creating a filter 
for filtering the spectra to emphasise higher resolu- 
tion parts of the unrequired components. 

19. A method according to claim 18 wherein the filter 
comprises a band pass filter. 

20. A method according to claim 18 or claim 19 wherein 
the filtering step also includes carrying out a least 
square fit of the filtered sample spectral data to the 
filtered unwanted component spectra. 

21. A method according to claim 20 including iterating 
the least square fit to remove those parts of the 
spectrum which have a poor fit. 

22. A method according to claim 20 or claim 21 includ- 
ing using the resulting co-efficients to compute an 
unfiltered unwanted component spectrum which is 
then subtracted from the measured sample spec- 
trum. 
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